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Abstract 

Isoscalar collective modes in a relativistic meson-nucleon system are investi- 
gated in the framework of the time-dependent Thomas-Fermi method. The 
energies of the collective modes are determined by solving consistently the 
dispersion relations and the boundary conditions. The energy weighted sum 
rule satisfied by the model allows the identification of the giant ressonances. 
The percentage of the energy weighted sum rule exhausted by the collective 
modes is in agreement with experimental data, but the energies come too 
high. 
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I. INTRODUCTION 



Renormalizable relativistic quantum field theories of liadronic degrees of freedom, called 
quantum hadrodynamics (QHD), have been studied for some time At the level of the 
mean-field theory (MFT) and one-loop approximation, these models have proven to be a 
powerful tool for describing the bulk properties of nuclear matter. The binding energy of nu- 
clear matter in MFT arises from a strong cancellation between repulsive vector and attractive 
scalar potentials. Such potentials are comparable to those suggested by Dirac phenomenol- 
ogy 113) Brueckner calculations 0], and finite-density QCD sum rules p. Therefore, it is 
not obvious that QHD would be able to reproduce the spectrum of finite nuclei, involving 
energies of the order of tens of MeV. However, it has been shown that it can realistically 
describe densities, single-particle energies and the spectrum of collective excitations of finite 
systems iiif- 

Collective modes of a relativistic many-body system are characterized as poles of the 
meson propagator. However, in the one-loop approximation, the meson propagators have 
also poles at space-like momenta, which arise from polarization effects of the Dirac sea 
p!0Hl3[ . While the existence of these poles does not rule out meson-nucleon field theories as 
useful descriptions of nuclear systems at low g, it may restrict the range of validity of several 
approximations to these theories. To avoid this problem, in this work we will study collective 
excitations of finite nuclear systems in a semiclassical approximation to the Walecka model. 

In refs. p!^,P^ a semiclassical approximation to the Walecka model was introduced to 



study collective modes in nuclear matter at zero and finite temperature. It was found that 
the results obtained are compatible with microscopic calculations of the meson propagators 
H'lll- want to generalize this semiclassical approach to the description of collective 



modes of finite nuclei by using a nuclear fluid-dynamical model which incorporates 

monopole and quadrupole distortions of the Fermi surface. This nuclear fluid-dynamical 
model has recently been applied with success to the description of temperature effects in 
collective excitations of finite nuclei [ pO[| . 

In ref. (which is a generalization of the works presented in refs. |0J^) isovector and 
isoscalar collective modes were calculated in the Walecka model, by introducing local hy- 
drodynamic variables to describe the nucleon fluids with the assumption of irrotational flow 
and in the limit of large masses for the vector mesons. As suggested in ref. , we lift these 
restrictions and in this work we calculate the isoscalar collective modes in the Walecka model 
in the framework of the time-dependent Thomas-Fermi method. 



In Sec. II we extend the formalism developed in refs. [p!8| , pJ| to the Walecka model. 
Collective modes are described by allowing the meson-fields and the nucleon densities to 
acquire a time dependence. The nucleon motion modifies the source terms in the meson 
field equations producing corresponding time-dependent changes in the meson fields. Since 
the nucleon dynamics is in turn specified by the meson fields, collective modes of nuclear 
motion arise naturally in this approach. In Sec. HI we derive the equations of motion, 
boundary conditions and orthogonality relations that the normal modes must satisfy. The 
dispersion relations, which solved consistently with the boundary conditions, determine the 
eigenvalues, are presented in Sec. IV. In this section the sum rule satisfied by the model is 
also given. We identify two rather collective monopole modes at 28 MeV and 35 MeV. These 
large values are expected since the isoscalar monopole excitation is a compression mode and, 
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therefore, its energy is related to the compressibihty of nuclear matter |2T|], which is known to 
be too low in the Walecka model. For the other multipolarities, we also observe that the most 
collective states come at higher energies than the experimentally observed giant resonances. 
It is true that our lowest modes coincide with the modes obtained by 0. However, these 
modes only carry a small percentage of the energy weighted sum rule and therefore should 
not be identified with the giant resonances. Finally, in Sec. V we give our numerical results 
and conclusions. 

II. FLUID-DYNAMICAL MODEL 

In a classical approximation to the Walecka model the energy of a nuclear system is given 
by H 



E = aI /(x, p, t) {[(p - g^Y)' + {M- g^a)']'/' + g^V^} 



{2nf 

+ \jd^x {III + Va ■ Va + m^a^) 

+ \jd\ [Ul^ - 2ny,a,K) + VV, ■ VV, - djVAVj + - V^)] , (2.1) 

where the distribution function, /(x, p, t), is restricted by the requirements 

^ = 4y^^/(x,P,t), (2.2) 

f (x,p,t)-/(x,p,t) = 0, (2.3) 
and its time evolution is described by the Vlasov equation 

f + {/,M = 0, (2.4) 



where h = y {p — gvY)"^ + (M — 5'sO")^+5'„Vo = e+gyVo is the classical one-body Hamiltonian 
and {,} denote the Poisson Brackets. 

The time evolution of the fields is given by 

VV + m^d = 5'sPs(x, t) , (2.5a) 



dt 
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d^Vc 



V72l 



-I 


dVo 


dt 


[ dt 
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(dVo 


dxi 


[ dt 



(2.5b) 



Qt2 - ^ + = 9vM^^ ^) + — I ^ + V ■ V ) , (2.5c) 



with 
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p,(x,t) = 4y ^^/(x,p,t) , (2.6) 

PB(x,t)=4|^/(x,p,t), (2.7) 



Using the Vlasov equation, Eg . (p^) , it can be easily shown that the four-current satisfies 



the continuity equation, and that the components of the vector field are related through |0 



d^V^ = . (2.9) 

Therefore, the second term on the right-hand side of Eqs. (|2.5b| , |2.5c| ) vanishes. 

In our calculations we will assume that the density of a spherical nucleus in the ground- 
state is constant inside the nucleus and zero outside, and is given by 

Po(r)=4/^/o(x,p). (2.10) 

with 

/o(x,p) = e[p^(r)-p2], (2.11) 

where Pri^) = Pf^['^~^]^ Pf is the nuclear matter Fermi momentum, and Rq is the nuclear 
radius. The ground-state distribution function /o is determined by the particle number A 
and by the minimization of the energy and the equilibrium nuclear matter density, pg, is 
calculated from equations ( |2.11| ) and ( |2.10| ) 

Po{r) =poQ[Ro-r]. 

Giant resonances manifest themselves as small amplitude highly collective modes. There- 
fore, they are described at the microscopic level by the RPA equations. In the classical limit, 
these equations are obtained by the linearization of the Vlasov equation. In this context we 
begin by expanding the distribution function around its equilibrium value /o(x, p): 

/(x, p, t) = /o(x, p) + {S, fo} + ^{S, {S, fo}} + ... , (2.12) 

where S'(x, p, t) is a generating function which describes small deviations from equilibrium. 

In its more general form, the distribution function, /(x, p, t), should include static as 
well as dynamic deformations of the nuclear system. For this reason we decompose the 
infinitesimal generator S'(x, p, t) into a time-even and a time-odd part 

^(x, p, t) = P(x, p, t) + Q(x, p, t) , (2.13a) 
Q(x,p,t) =g(x,-p,t) , (2.13b) 
P(x,p,t) = -P(x,-p,t). (2.13c) 
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The time-even generator, Q(x, p,t), takes into account the dynamic deformations. The 
static deformations are described by the time-even distribution function, which includes the 
fields responsible for the deformations of the Fermi surface. In the present approach, it is 
expressed in terms of the time-odd generator P(x, p, t) 

/^(x, p, t) = /o(x, p) + {P, /o} + ^{P, {P, /o}} + ... 

= e[A - /io(x, p) - W{^, t) - ]^PiPjXiA^, t)] . (2.14) 

The scalar field, W{'x, t), is related to the deformations which preserve the spherical form of 
the Fermi surface. The tensor field, Xjj(x, t), introduces deformations in the Fermi sphere. 
Hopefully, the scalar and tensor fields will provide an adequate description of the monopole 



and quadrupode deformations of the Fermi sphere. In Eq.( |2.14| ), /?.o(x, p) = yP^ + M*^(x)-|- 
^^V2^(x), with M*(x) = M — cro(x), and cro(x) and lo°(^) ^^^^ respectively, the equilibrium 
values of the fields a and Vq. The Fermi momentum is related to A through 



A = iplir) + M*2(x) + gXi^) =eF + gXi^) ■ (2-15) 
The introduction of the generator Q(x, p, t) destroys the time reflexion invariance of the 



equilibrium distribution function. It will allow for the appearance of transverse flow 
the nucleus. The simplest choice which includes this possibility is given by |T 



m 



Q(x, p, t) = ?/^(x, t) + ^piPj(f)ij{:>c, t) , (2.16) 

where ilj{:K,t) and (l)ij{:s.,t) are, respectively, scalar and symmetrical tensor fields. 

The time evolution of the generator S and the field fluctuations are determined by the 
appropriate Lagrangian. For small deviations from equilibrium it is enough to consider the 
quadratic Lagrangian 

L^'^ = 2 / fo{S, S} + J d'xU^a + J d'xUyV, - E^'^ . (2.17) 

Using the ansatz Eqs. (|2.12|) , (|2.14 ) and (|2.16|) , decomposing the boson fields into a static 
(ground-state) contribution and a small time-dependent increment and imposing the barion 
number conservation, we get 



(27r 



, PfPo ( Sij \ ( ] 5ij ] 
H Yq' \ '^''■^ ^Xfcfe I I (Pij -^(Pkk 



+ /rfS.5RpoU + ^0j , (2.18) 
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£;(2) 



(fx 



■Z V 



Po 
2eF 



20 I 2 



2 4 



280e, 



+ 2V0ij ■ V0ij + 4:di(j)ijdk(l)kj + 4:dk(l)ijdj(f)ik)] + J d^x [{gsM*6a 



+ ^Jd^x [Ul + V6a ■ V6(x + {ml + Am^)((5a)2] + ^ J d^x [H' 
- 2UvASVo + V5y. ■ V6Vi. - d,6VASVj + {ml + Aml){6Vf 
-ml{6Vo)^] + J d'E.6R{gsPso^a-g^PoSVo) . 



(2.19) 



The surface integrals in the above equations take into account possible surface displacements 
parametrized by a vector field, 5R(x). Our choice of the even distribution function allows 
exphcitly for this effect. In Eq. (|2l^) , Ami = 9s§m^ and Ami = 91po/^f- 



III. EQUATIONS OF MOTION, BOUNDARY CONDITIONS AND 
ORTHOGONALITY RELATIONS 



The equations of motion and boundary conditions that specify the dynamics of the fields 
are obtained from Eq. ( |2.17|) through the Euler-Lagrange equations. We get 



6& = n. 



(3.1a) 



- V'6a + {ml + Aml)6a = -g,M* (^W + ^x^^ , , 



(3.1b) 



6V, = Uv, - d,6Vo 



Pi 



6V, - V'6V, + {ml + Aml)6V, = -^po + ^^^'t>n + 25j0u) > 
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(3.1c) 
(3.1d) 



(3.1e) 



(3.1f) 
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(3.1g) 



73 



2 \ 2 2 2 



+ d,6Vi) + 
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Pf 
35el 



^ij ( |V^0A:fc + dkdi(f)ki ) + VVij + ^i^j^fcfc + '2didk4>kj + 2djdk^Pki 



(3.1i) 



It is worth mentioning that Eqs. (|3.1a|) to (|3.1i| ) are vahd only in the interior of the nucleus. 
Therefore, we replace pi?, ep and po in these equations by their equilibrium values. At the 
surface, the variational fields satisfy the following boundary conditions 



Xk{dk5a + gsPso^Rk)\r=Ro = 



(3.2a) 



Xk{.dk5Vi - di5Vk%^^^ = 



(3.2b) 



XkidkSVo + 5Vk + g^poSKk) 



r=Ro 



(3.2c) 



r=Ro 



(3.2d) 



Xk 



' • 1 p2 

SRkSij + —{dki^Sij + diijjSjk + djipSik) + ^{di(j)ki5ij 



+ dj(l)iiSik + di(j)ij6ik + di(j)ii6jk + -dk(f)ii5ij + dk(j)ij + dj(j)ik + di(j)jk) 



+ ^{SVkSij + SViSjk + SVjSik) - i^iSjk + ijkk] 



r=Ro 



(3.2e) 



10 Po 



r=Ro 



(3.2f) 



In order to ensure that the current density is not singular at the surface, the following 
boundary condition has also to be imposed [ITU 



Xk(pkj\ 



r=Ro 



. 



(3.3) 
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In Eq. (|3.2eO , ^ is a vector Lagrange multiplier that takes into account the restriction 

We look for normal-mode solutions where all the fields oscillate harmonically in 
time. This means that the fields are described by a superposition of the real parts 
of {n^"), V'^"^ ^V("), i5a^'^\ i5Vi''\ iW^'^\ i(5R("), ill.^y\ ixjj^} exp"*^-*, where all the 
quantities within the braces are only functions of x. This normal-mode analysis leads to the 
RPA coupled equations for the eigenmodes: 

uja^^^ = n(") , (3.4a) 
- ul5a^-^ - V^5a(") + {ml + Am^)5aW = -g^M* {^^W^^^ + |x!r^) , (3.4b) 

- ujjvt^ = ni,^) - d,6Vt^ , (3.4c) 
- ul5V,^-^ - VHvt^ + {ml + /\ml)vt^ = -^p, (d,^^^^ 

+ |M;^+25.4"0) ' (3.4d) 
- ul5V,^-^ - VHvt^ + mlSV^-^ = -epg. {^W^^^ + , (3.4e) 



n ( W^-^ + ^xlr^) = ^V^^(") + || (W.? + 2d,d,4f) + ^-^dM""^ , (3.4f ) 



)- 


















- un^f = ^ {i ^ j) , (3.4h) 



ep I ep 



It is clear form Eq.( |3.4h| ) that Xij and (pij are canonically conjugate fields. 

The solutions of the the above equations satisfy the following orthogonality relation 
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IV. DISPERSION RELATIONS AND SUM RULES 



A. Dispersion Relations 



The eletric modes are described by the same kind of solutions as constructed in ref. [jlQ 
i. e., by two kinds of transverse fields 



= - 5,,V^)f - [9,(V X 1), + d,{V X I),] 



- [(V X I),(V X I), + (V X I),(V X \),]}ji{hr)Yio , (4.1) 



[<P^J]2 = X 1).- + X Wiik2r)Yio , (4.2) 

and by one longitudinal tensor field 

[<P^J]s = (d,d, - I^V^^ Ji{hr)Yio . (4.3) 

The advantage of using the above combination of the four linearly independent angular 
tensor functions: didjVio, SijYiQ, {xidj+Xjdi)Yio, and XiXjYiQ, is that all solutions given above 
are traceless. In particular, the transverse fields also verify the relations 

di[<Pi,]^ = and did, [<Pij]^ = . (4.4) 

For each multipolarity, all scalar fields are proportional to ji{kr)Yio, and the vector fields are 
combinations of two linearly independent vector functions: di{ji{kr)YiQ) and (Vxl)jj;(A;r)y;o- 
Using these combinations in Eqs. (|3.4a|) to ( p.4i|) it is straightforward to show that the 
transverse solutions do not couple to the scalar fields, and one has [5cr]i,2 = [5Vo]i,2 = 
[W^]i,2 = [V']i,2 = [no-]i,2 = 0. For solutions of kind 1, the vector fields are also zero: 
[SVi\i = [IIvJi = [SRi]i = and the dispersion relation for this particular solution is given 
by 

ul = ^^kl . (4.5) 



This is the same relation as obtained in ref. [jT9[. This should be expected since the meson 



fields, which are the new ingredients in the model used here, do not couple to the solution 
of kind 1. 

For solutions of kind 2, we still have [SRi]2 = 0, since, from Eq. (|3.2a|) , the vector field 
5R is directly related to the scalar field Sa. However, the vector fields [6Vi\2 and [nyj2 are 
coupled to the tensor fields. We get 

mh = ^^fA [0^2 = -^^Vlkliy X \U{k2r)Y,, , (4.6) 

where 
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and m*^ = + Am^. Using the solutions of kind 2 in the normal mode equations we get 
(-^ - ^) U - kl - mf) = , (4.8) 

which give us two different solutions for For gy = 0, one of the solutions is exactly 



the same which is obtained in ref. |T^. This solution is now modified and a new solution 
appears, due to the coupling between the vector meson field and the fields introduced to 
describe the nuclear deformations. 

The longitudinal solutions, [4>ij]^, couple to all other fields and give 

li'h = f(h)ji(hr)Ym , (4.9a) 



\W], = , (4.9b) 



= m*^) ^^^' = ^(^3h(A^3r)y,o , (4.9c) 

^ ['^n ~ ~ ) 



= .^J'^'p\^) ^^^' = ^o(^3)jKA:3r)l^o , (4.9d) 



= ^ ( fih) - dMhr)Y,o) , (4.9e) 



9v V 15 
Xidi[6a] 



Xi[6R, 



i\3 



9sPsO 



(4.9f) 



r=Ro 



plus the corresponding solutions to the canonically conjugated fields. In the above equations 
we have introduced the functions 

r.j. 2p%k\l + Gy{k))Gos{k) 

> 15{3elujl-plk^{l + Gy{k))Gos{k)) ' ^ ' ' 

and 

G^jk. ^ 1 _ ( alel \ ^ 

-n'^ep \ijJn — k"^ — rnl — k"^ — m*^ ' 



with Gy{k) defined in Eq.(p|). 

The dispersion relation obeyed by these solutions is 
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9 



T 



[9yA< - kl - m-) - giM*\ul - - ml)] . 



(4.12) 



There are four solutions of kind 3, two more than the number of this kind of solutions found 
in [|^. This should be expected since, besides the vector meson field, the scalar meson field 
also couples to the longitudinal solution [(f>ij]3- It is easy to show that for Qs = and Qv = 



one recovers the two solutions of ref. 19 



Therefore, the Walecka model leads to the appearance of 7 different values for k for a 
fixed frequency u, in contrast with the model of ref. [|l9l, which gives only 4 different values. 
There is still a fourth kind of solution for the tensor fields, which can be chosen to be 



coupled to the scalar fields 



and to the meson fields 



SijF{r)Yio 



Pf 



[6aU = [SVoU = [SV^4 = , 



(4.13) 



(4.14) 



(4.15) 



where F{r) is an arbitrary function. This solution is not trivial because of the boundary 
condition Eq.(|3.3|). 

The general solution, for each normal mode, is a linear combination of the eight particular 
solutions: 



= Cl[(j)ij{kir)]i + J2 C2n[<Pij{k2nr)]2 + C3„[0ij(A;3„r)]3 + C4[0y(r)]4 



(4.16) 



n=l n=l 

with similar expressions for the other fields. 

To avoid zero frequency modes linked to the surface motion, we introduce in the model 
a surface energy which, in a classical approximation, is given by 



^S = ^(/(^ + l)-2)/rfS.5R5R.n 



(4.17) 



where a sup is the surface tension coefficient. This term does not alter the equations of motion 
and, therefore, the dispersion relations. It only changes the boundary condition Eq. (|3.2fl) to 



10 po ^Rm 



2)n ■ (5R 



. 



(4.18) 



r=Ro 



Using the general solutions in the boundary conditions Eqs. (|3.2a|) to (|3.2e| ), Eq. (|3.3| ) and 
Eq.( [4.18D we get the equations ( |6.1a| ) to ( |6.1h| ) given in the appendix. The eigenvalues are 
determined by solving consistently the dispersion relation equations, Eqs. ([4.5|) , (f4.8|) and 
( [4.12| ), subjected to the boundary conditions. 
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B. Sum Rules 



Sum rules can be regarded as a test to the validity of a particular nuclear model. Suppose 
that a nucleus is excited from its ground state |0) to an excited state |n), with an energy 
En, due to interactions with an external field. One can define momenta, weighted in energy, 
of the excitation strength distribution 



ruk 



Y.{Er,-Eor\{n\dm 



(4.19) 



where O is the one-body hermitian operator, responsible for the excitation. In the above 
expression, k = 0,±1,±2,... and \n) stands for a set of eigenstates of the hamiltonian of 
the system. A sum rule is obtained when it is possible to relate a momentum with a known 
quantity. 

The energy weighted sum rule (EWSR), mi, is obtained through the calculation of the 
expectation value of a double commutator 



m, = Y.{En - E,)\{n\dm' = -(0|[6, [if,6]]|0) . 



(4.20) 



In the present problem, the general solution for the variational fields is given by the real 
part of 

/ iiy(")(x) \ 
i5cT(")(x) 



^(x,t) = E' 



-iuJnt 



(4.21) 



i5R(")(x) 
zn(")(x) 

^Xii (x) 

0;;Hx) 

?/^w(x) 

nW(x) 

V (5V(")(x) ) 

where the coefficients a„ are determined by the initial conditions. In order to derive the 
EWSR for the electric modes we consider the following initial condition 



V^(x,0)=D(x) 



(4.22a) 



%(x, 0) = x^i(x, 0) = W^(x, 0) = 5R(x, 0) = 5V(x, 0) 

= (5a(x, 0) = 5Vo(x, 0) = n,(x, 0) = Hy^ (x, 0) = 



(4.22b) 



with -D(x) to be specified. We then expand the fields '?/'(x, 0), no-(x, 0), 0jj(x, 0) and 5V(x, 0) 
as V9(x, 0) = Z^n '^nV^*'"''(x), whcrc, from the orthogonality relation Eq.(|3.5|), we get 



W^-^ + ^xSr^) D{^)+p,J dS.5R(")D(x) 



(4.23) 
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The coefficients are related to the expectation value of the transition operator, a„ = 
v/2(n|6|0). Therefore, the EWSR can be written as 

mi = Y.\an\^uOn = 2E^^\ (4.24) 

n 

and, for the initial condition given in Eq.(4.24), the EWSR reads 

\an?0Jn = f d^x^VD.VD . (4.25) 

V. NUMERICAL RESULTS 

We have performed our calculations with two different sets of the mean-field values of 
the parameters in the Walecka model: 

I. gl = 122.88,^2 ^ 169.49, = 1.3fm-\ M*/M = 0.522 

II. c/2 = 91.64, c/2 = 136.20, = lA2fm-\ M*/M = 0.556 

where M = 938 MeV and M* is the effective mass. The effective mass and the Fermi 
momentum indicated for each set correspond to the values at which saturation of nuclear 
matter is obtained with an energy per nucleon E/N = —15.75 MeV, using rris = 550 MeV 



and = 783 MeV. The surface tension, from the liquid drop model is asup = 1-017 
MeV/fm^. The results were calculated for a nucleus with A = 208. The radius Rq is obtained 
from the value of pp corresponding to the chosen set of parameters. 



For the excitation operator introduced in |4.22a| we will use 

D(x) = rVoo, / = , (5.1a) 

= r'Yio, l>2. (5.1b) 

Table I shows the energies of the normal modes together with the corresponding per- 
centage of the exhausted energy weighted sum rule (EWSR), for the two sets given above 
and for different multipolarities. The EWSR is fragmented over the whole range of ener- 
gies and only the nuclear modes which exhaust more than 0.1% of the sum rule are given. 
The distribution of the EWSR between the nuclear modes and the mesonic modes (energies 
larger than the meson masses) agrees with the results obtained in []T1[, where it is shown 



that in infinite nuclear matter and for small momentum transfer about 62% of the EWSR 
is exhausted by the continuum nuclear modes and about 38% by the vector meson modes. 
For instance, for I = 2+ and for set I, we find a vector meson mode at hwi = 984.56MeV 
which exhausts 27.30% of the EWSR. The other mesonic modes are not as collective as 
this one and are distributed over a large range of energies. This pattern is reproduced for 
the two sets of parameters and for all multipolarities. The EWSR is fulffiled considering 
all the nuclear and mesonic modes. In non-relativistic calculations using the same nuclear 
fluid-dynamical model used here |ll8|-pO[], the mesonic modes are not present and, therefore. 



the EWSR is distributed only through the nuclear modes. From this Table we can see that 
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for set II the collective modes come at a slightly higher energy than in set I and that the 
strength is more concentrated at higher energies. 

In table II we give for set I and for l'^ = O"*", 2"*", 3~, 4"*" the energy of the normal modes 
with energy below 100 MeV (first column) and the corresponding percentage of the energy 
weighted sum rule (EWSR) (third column). In the second column we present a renormalized 
percentage of the EWSR, renormalizing the strength distributed among states with energy 
below 130 MeV to 1. The renormalizing factor is m'i{l) = 0.56, 0.60, 0.61, 0.56, respectively, 
for / = 0, 2, 3, 4. This is done so that we can compare more easily the results obtained in the 
present work with previous results obtained in a non-relativistic fluid-dynamical model, |]20| , 
(columns 4 and 5) and experimental data (columns 6 and 7) Looking at the modes with 
energy below 100 MeV, we may immediately conclude that there is a certain correspondence 
between the states obtained in the present approach and the ones of 1^^, if we identify the 



states by the percentage of the exhausted EWSR. However, the corresponding states come, 
in the present relativistic approach, at higher energies. For instance, the quadrupole low 
lying mode and giant resonance come, respectively, at 10 and 20 MeV and exaust 8% and 
77% of the EWSR while the experimental modes come at 4 and 11 MeV and exhaust 15% 
and 70% of the EWSR. Another possible way of identifying the modes is done by comparing 
the current transition density (2.8) and the transition density (2.7) for these two modes with 



the ones of ref. In figure 1 and 2 we plot j+, j_, jdiv (arbitrary units) defined by the 
equations: 

j(r) = j+(r)YM+i,o(fi) + j-(r)Y,,,_i,o(f^), 
V • j(r) = jdivir)Yio. 

The function jdiv is related to the transition density 6p (p = —V ■ j). For the 10.03 MeV 
mode, j+ and j_ have oposite signs and jdiv is close to zero, characteristic of a surface mode. 
These are typical properties of a low lying mode. For the 20.15 MeV mode, j+ and j_ have 
the same sign and jdiv comes diferent from zero for t/Rq > 0.5. This behaviour is closer 
to the behaviour expected from a giant resonance. We conclude the identification we have 
done is correct. 

We note that our modes with the lowest energy have energies similar to the ones obtained 
in ref. , however, these are not the states that exhaust the largest percentage of the EWSR 
and, therefore, they should not be identified with the giant resonances. The breathing mode 
comes at a very high energy, but this was expected owing to the high incompressibility of 
the model. 

While in ref. only the lowest modes were determined, we have found all the modes 
that exhaust a significant fraction of the corresponding EWSR (which we also derived). 
Furthermore, we have shown that the lowest modes are not the most collective ones. 

From the present results we conclude that the dynamical properties of the nuclei are 
not so well described by the Walecka model as the static properties such as densities and 
single particle energies. In our calculation we have taken for the ground-state of the nucleus 
a Slater determinant derived from a square-well instead of the the self-consistent ground- 
state. We believe, however, that for large nuclei such as the ^osp]^ nucleus this is a good 
approximation which allows us to obtain analytical expressions for the equations of motion 
and the boundary conditions. 
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VI. APPENDIX 



Using the general solutions in the boundary conditions Eqs.( |3.2aD to ( |3.2e| ), Eq.(|3.3|) and 
Eq.( [4.18| ) we get the following equations: 



n=l 



(2 - /(/ + l)){rdr + l)ciji{kir) + ^ [2l{l + 1) - 2rdr -2-k 

4 

+ ^(rdr - l)c3„ji(A;3„r) 



2 2 
2n^ 



C2njl{k2nr) 



n=l 



r=Ro 



(6.1a) 



1(1 + 1)(2 - /(/ + l))cijiihr) + 2l{l + l){rdr - I)c2njiik2nr) 

2 



n=l 



+ E 

n=l 



/(/ + 1) _ 2rdr - ~-kiy 



csnjiihnr) + r CiF{r) 







r=Ro 



(6.1b) 



ky{rdr - 1) + Q{rdr + 1) - ^-!^-^rdr 



c^ji{k^r)~2Y,[&{rdr + l) 



n=l 



+2kiy - 3l{l + 1)1 C2nJi{k2nr) + 3j2{rdr - 2)CsnJl{k3nr) 



n=l 



r=Ro 



(6.1c) 



/(/ + 1) [ky{rdr - 1) + 3rdr + 12- + 1)J ciji{kir) 

+ ^ /(/ + 1) (hiy + Qrdr - 12) C2nJl{k2nr) 
n=l ^ 



+ E [^^L^'-^, + 3 (-3ra, + /(/ + 1) - A;: 



n=l 



2 2 



c-in3i{kznr] 



r=Ro 



(6.1d) 



2 

E 

n=l 



Gv{k2n)k2nr C2nji(^2n'^) 



(6.1e) 



E ^"^;'^^^-) p^4/(/ + I)c2aik2.r) + E [K,(A:3.) 

UJnGv{kzi) 



1^2 7,2 



fik 



3i) 



15 









rdrCsijiik^ir) 


/ 9sPsO 




15 





r=Ro 



(6.1f) 



2 ^2 l2 4 



15 



9sPsO 



rdrC3iji{k3ir] 







r=Ro 



(6.1g) 



Po 



a 



sup 



■(/(/ + l)-2)a(A;3,)9, 



r=Ro 



(6.1h) 



with the functions a{k) and Vo(/i;) defined in Eqs. ([4.9c| ) and (|4.9d|) . 
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TABLES 



TABLE I. Energies and fractions of the energy-weighted sum rule for different multipolarities 
and different sets of parameters. 

TABLE II. Comparison between the energies and fractions of the energy-weighted sum rule 
obtained in the present work (first column), in [20] (second column) and experimental data (third 
column) [24]. 
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FIGURES 

FIG. 1. j+ (full- line), j_ (dashed-line) and jdiv (dash-dotted-line) in arbitrary units for the 
= 2+ E=10.03 MeV mode. 

FIG. 2. j+ (full-line), j_ (dashed-line) and jdiv (dash-dotted-line) in arbitrary units for the 
= 2+ E=20.15 MeV mode. 
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Table I 



7 


1 


11 








/icc;i(MeV) 


mi(%) 




28.56 


14.06 


37.27 


8.55 


z 


35.50 


27.46 


46.12 


30.34 




50.95 


2.30 


61.59 


2.13 




68.25 


1.32 


81.14 


0.25 





71.29 


5.94 


86.73 


8.59 




88.51 


0.58 


98.92 


0.71 


0^ 


105.11 


1.07 


133.51 


3.22 


o 


107.87 


2.57 






2^ 


10.03 


4.67 


11.90 


1.46 


2^ 


20.15 


45.32 


28.07 


42.51 


2^ 


28.32 


0.78 


33.77 


5.12 


2+ 


35.32 


3.95 


39.94 


0.48 


2^ 


35.82 


0.31 


42.33 


4.64 


2^ 


49.35 


0.52 


59.72 


0.35 


2+ 

i 


64.20 


1.93 


75.59 


0.96 


9+ 

^8 


69.91 


0.10 






37 

1 


12.93 


11.42 


15.86 


0.32 


32 


14.44 


1.09 


17.77 


5.71 


3s 


32.84 


35.48 


42.14 


35.24 


37 


37.20 


0.70 


44.95 


4.09 


37 


42.33 


0.08 


47.65 


2.17 


36 


45.03 


3.90 


51.98 


7.00 


3f 


57.55 


0.64 


68.89 


0.53 


4^ 


18.06 


13.41 


22.35 


0.69 


4^ 


20.46 


1.37 


24.93 


7.21 


4^ 


43.94 


16.44 


54.32 


10.32 


41 


45.91 


14.54 


57.52 


20.27 


4^ 


48.84 


2.05 


62.47 


11.42 


4^ 


54.83 


3.47 


73.41 


0.11 


4^ 


65.46 


0.68 


77.68 


0.64 
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Table II 



experimental p4 



present 



2C 



huJiiMeV) mi{%)/m\{l) 



huiiMeV) mi{%) 



^2 

ot 

% 

ot 

total 



2t 
2t 
2t 



9+ 
9+ 
total 



28.56 
35.50 
50.95 
68.25 
71.29 
88.51 



25.23 

49.28 
4.14 
2.37 
10.66 
1.03 
99.99 



14.06 
27.46 
2.30 
1.32 
5.94 
0.58 
55.72 



15.87 

18.95 
28.14 
36.83 
41.29 



95.15 
2.26 
0.03 
0.03 
1.46 

98.98 



13.9 



100. 



10.03 
20.15 
28.32 
35.32 
35.82 
49.35 
64.20 
69.91 
87.43 



7.67 
76.99 
1.29 
6.35 
0.53 
0.88 
3.15 
0.17 
0.29 
99.88 



4.56 
45.79 
0.77 
3.78 
0.31 
0.52 
1.87 
0.10 
0.17 
59.46 



3.73 

11.70 

17.45 

20.54 

21.12 

27.30 



30.90 
64.19 
2.17 
1.10 
1.00 
0.06 



99.32 



4.09 
10.9 ± 0.3 



15. 
70.0 



3r 

33 
34 

35" 

36 

3f 

38- 

39 
total 



12.93 
14.44 
32.84 
37.20 
45.03 
57.55 
78.43 
82.04 
95.87 



18.74 
1.71 

58.70 
1.15 
6.06 
1.05 
2.96 
3.10 
0.46 

100.00 



11.43 
1.04 
35.79 
0.70 
3.69 
0.64 
1.81 
1.89 
0.28 
61.00 



2.92 

8.43 

18.53 

22.80 

26.87 



34.10 
0.29 
43.44 
10.88 
5.18 



2.61 

18.4 ± 0. 
21.8 ± 0. 



33. 

36. 
27. 



97.64 



4r 

4^ 



total 



18.06 
20.46 
43.94 
45.91 
48.84 
54.83 
65.46 



23.10 
2.35 

29.21 

25.77 
3.62 
5.81 
1.20 

100.00 



13.06 
1.33 
16.51 
14.57 
2.05 
3.28 
0.68 
56.52 



4.51 

12.26 

23.36 

27.64 

29.67 

33.45 

35.38 



34.10 
2.05 
22.39 
8.86 
17.10 
8.45 
4.18 
97.13 



4.32 
12.0± 0.3 



10 ± 3 
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